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We  theoretically  examine  the  idea  of  using  the  Electrohydrodynamic  (EHD)  ionic  wind 
pump  effect  in  thruster  applications.  This  idea  has  been  under  discussion  for  many  decades. 
However,  there  is  still  no  commonly  accepted  view  on  EHD-thrusters,  e.g.,  whether  EHD- 
thruster  is  worthwhile  or  not,  what  level  of  thrust  and  thrust  efficiency  can  be  obtained  from 
EHD-thrusters,  etc.  In  this  paper  we  present  a  simple  one-dimensional  model  of  an  ideal 
EHD-thruster  for  calculating  thrust  efficiency  and  thrust  of  EHD-thrusters.  We  also 
calculate  the  maximum  current  that  can  be  achieved  for  an  ideal  EHD-thruster  at  a  given 
voltage.  This  allows  us  to  calculate  the  maximum  thrust  that  can  be  obtained  from  the 
thruster  and  the  corresponding  thrust  efficiency.  We  show  that  with  an  increase  in  the 
voltage,  the  maximum  trust  and  the  corresponding  thrust  efficiency  move  in  opposite 
directions:  the  thrust  efficiency  decreases,  while  the  thrust  increases.  Our  model  illuminates 
the  physical  limitations  of  EHD-thrusters  and  provides  reasonable  estimates  of  the 
performance  of  real  EHD-thrusters. 


Nomenclature 


B  =  calculated  constants 

=  magnitude  of  electrical  field 
i  =  magnitude  of  electrical  field  at  the  anode 
e  =  electron  charge 

H  =  altitude  above  sea  level 

i  =  index  corresponding  to  z-th  spatial  step 

j  =  electrical  current  density 

Jmax  =  maximum  space-charge-limited  current 

L  =  the  discharge  length  of  the  thruster 

ma  =  mass  of  neutral  gas  molecule 

nij  =  mass  of  an  ion 

na  =  neutral  gas  number  density 

na  o  =  neutral  gas  number  density  at  the  anode 

«,■  =  ion  number  density 

P  =  gas  pressure 

T  =  thrust  per  unit  area  of  the  discharge 
Ta  =  temperature  of  neutral  gas  flow 

Tj  =  ion  temperature 
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rmax  =  maximum  thrust  per  unit  area  of  the  discharge 
U  =  discharge  voltage 

Uq  =  a  given  discharge  voltage 

u  =  velocity  of  neutral  gas  flow 

Uq  =  incoming  velocity  of  neutral  gas  flow  at  the  anode 

M|  =  exit  velocity  of  neutral  gas  flow  at  the  cathode 

V drift  =  i°n  drift  velocity 

V drift  o  =  t°n  drift  velocity  at  the  anode 
Vja  =  thermal  velocity  of  neutral  gas  flow 

Wj.  =  gain  of  kinetic  energy  of  gas  flow  per  second  per  unit  area  of  the  discharge 
Wj  =  electrical  power  deposited  in  the  discharge  per  unit  area  of  the  discharge 
v  =  axis  directed  along  the  thruster 

A Ta  =  increase  the  gas  temperature  in  the  thruster  due  to  Ohmic  heating 

AV  =  delta-V 

Aa  =  spatial  step 

£q  =  permittivity  of  free  space 

lw  =  ion-neutral  collision  length 

Mj  =  Hi-na 

/jj  =  ion  mobility 

tui  =  ion-neutral  collision  time 

X  =  thruster  efficiency 


I.  Introduction 

Electrohydrodynamics  (EHD),  or  alternatively  electro-fluid-dynamics  (EFD)  or  electrokinetics,  refers 
to  motion  of  electrically  charged  fluids  under  the  influence  of  applied  electric  fields.  In  EHD-thrusters  the 
applied  electric  field  accelerates  the  electrically  charged  particles  while  collisional  coupling  between  the 
charged  particles  and  neutral  particles  is  used  to  convert  the  directed  energy  of  the  charged  particles  into 
neutral  gas/fluid  particle  flow  without  the  need  for  moving  parts.  Fig.  1.  Since  the  momentum  transfer 
between  ions  and  neutral  particles  is  much  more  efficient  than  between  electrons  and  neutral  particles, 
neutral  gas  flows  produced  in  EHD-de vices  are  commonly  referred  to  as  ion  wind  (also  ionic  wind,  ion- 
driven  wind,  corona  wind,  electric  wind,  and  others)  flows.  The  EHD  ionic  wind  pump  effect  was 
discovered  1709  by  Francis  Hauksbee  [1]  and  later  investigated  by  numerous  scientists  including  Newton, 
Faraday,  and  Maxwell;  a  history  of  EHD  ionic  wind  pump  effect  can  be  found  in  excellent  review  of 
Robinson  [2],  The  predominantly  neutral  gas  flows  produced  using  the  EHD  pump  effect  have  been 
successfully  incorporated  into  a  variety  of  practical  devices.  The  most  prominent  application  is  currently 
ionic  air  purifiers  such  as  the  Ionic  Breeze  [3],  Such  devices  use  the  ion  wind  effect  to  drive  untreated  air 
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through  the  device  and  hit  a  collection  plate  where  the  unwanted  particles  stick  to  its  surface.  Another, 
more  recent,  application  is  EHD  laptop  cooling  fans.  Jewell-Larsen  with  his  colleagues  [4]  have 
demonstrated  that  such  devices  can  achieve  similar  cooling  performance  as  the  existing  high  RPM  cooling 
fans,  but  at  reduced  noise  levels.  In  1961  Steutzer  [5]  investigated  theoretically  and  experimentally 
operating  the  EHD  ion  wind  devices  backwards  as  a  high  voltage  generator,  where  the  ion  current  is 
generated  by  friction  with  liquid  flow. 

EHD  based  propulsion  systems  were  proposed  and  investigated  in  the  early  1960s.  The  most 
prominent  effort,  called  the  Ionocraft  and  led  by  Major  de  Seversky,  made  the  front  cover  of  Popular 
Mechanics  in  August  1964  [6],  The  article  described  the  technology  as  a  "magic  carpet"  that  was  viable 
over  the  entire  range  of  speeds  up  to  and  including  jet  aircraft  speeds.  It  also  claimed  that  it  would  become 
more  efficient  than  propellers  or  jets  for  aircraft  propulsion  and  be  able  to  operate  at  altitudes  as  high  as  60 
miles.  One  envisioned  application  was  an  Ionocraft  as  large  as  a  city  block  used  as  an  antimissile  platform 
that  would  indefinitely  stay  in  the  upper  atmosphere.  At  the  time  of  the  article,  however,  the  engineers  [6] 
had  built  only  a  rudimentary  device  out  of  balsa  wood,  aluminum  wire,  and  an  electrical  tether  down  to  a 
high-voltage  power  supply.  They  also  achieved  very  limited  performance:  it  took  90W  (30kV,  3mA)  to  lift 
2oz.  At  that  time  they  stated  that  their  active  work  was  focused  on  increasing  the  energy  efficiency. 

The  serious  theoretical  and  experimental  investigation  of  EHD-thrusters  began  with  the  work  of 
Christenson  and  Moller  [7]  in  which  they  developed  a  ID  model  of  EHD  propulsion  systems  and  compared 
it  with  experiments.  Experimentally  they  achieved  a  thruster  efficiency  (conversion  efficiency  from 
electrical  energy  to  fluid  kinetic  energy)  that  was  about  one  half  of  model  predicted  values.  They  concluded 
that  approximately  90%  of  the  input  energy  ended  up  as  flow  thermal  energy  and  that  the  ion  mobility  of 
air  would  have  to  be  reduced  by  about  2  orders  of  magnitude  for  realistic  applications.  In  their  theoretical 
analysis  the  authors  have  assumed  that  the  ion  mobility  is  constant  and  independent  of  gas  density. 
Although  in  general,  this  assumption  is  inaccurate  (it  works  only  when  AH  — >0),  their  model  gives 
reasonable  estimates  of  the  magnitude  of  the  thrust  and  the  trust  efficiency.  We  also  would  like  point  out 
on  work  of  Robinson  [8]  (which  was  closed  to  the  public  for  a  long  time),  where  he  had  already  in  1961 
developed  a  very  similar  model  as  Christenson  and  Moller  [7]  for  EHD-devices  in  which  he  has  also 
included  losses  specific  to  the  corona  discharge.  Robinson  came  to  the  same  conclusion  as  Christenson  and 
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Moller  [7]  that  the  conversion  efficiency  of  EHD  ion  wind  blowers  (or  EHD-thrusters)  is  very  low.  He  also 
has  pointed  out  that  in  certain  gases  the  generation  of  extraneous  gases  in  the  corona  discharge,  e.g.,  ozone 
in  air,  can  be  a  significant  limitation  of  using  corona  discharge  in  EHD-blowers.  In  1986  Bondar  and 
Bastien  [9]  analytically  estimated  the  thruster  efficiency  and  showed  that  an  increase  of  fluid  velocity  leads 
to  a  higher  x  ■  They  have  experimentally  achieved  a  conversion  efficiency  of  up  to  7.5%  by  operating  with 
an  incoming  flow  velocity  of  50  m/s,  which  was  close  to  their  estimated  value.  In  2005  Singhal  and 
Garimella  [10]  constructed  a  transient  3D  coupled  EHD/Navier-Stokes  code  and  used  it  to  model  the  effect 
of  incoming  velocity  of  neutral  gas  on  efficiency  of  EHD-pumps  (or  EHD-thrusters).  As  Bondar  and 
Bastien,  they  concluded  that  significant  increases  in  conversion  efficiency  of  EHD-pumps  are  possible  by 
increasing  the  incoming  velocity  into  the  EHD-pumps.  In  recent  work  [11]  the  authors  have  shown  that  the 
exhaust  velocity  can  be  increased  by  adding  multiple  stages  of  EHD-thrusters  in  series.  They  demonstrated 
a  maximum  axial  velocity  for  a  single  stage  2.5  cm  diameter  device  of  4.5  m/s.  They  also  tested  up  to  7 
stages  in  series  yielding  a  maximum  velocity  of  over  7m/s,  but  at  a  conversion  efficiency  of  around  0.1%. 

Thus,  the  primary  limitation  of  EHD  based  propulsion  systems  is  that  they  are  typically  very  inefficient 
in  terms  of  thrust  efficiency  when  operating  on  gases,  roughly  1-10%,  while  competing  technologies  like 
fans  can  achieve  %  of  60%-70%  [8],  There  are,  however,  other  factors  that  could  limit  the  applicability  of 
EHD  based  propulsion  systems.  One  of  them  is  that  the  ion  current  is  charge  limited  [7,  8]  that  leads  to 
small  ionization  fractions  in  EHD  thrusters,  typically  in  the  order  of  10'10  at  one  Atm  air  pressure.  This 
indicates  that  space  charge  built  in  the  EHD  thruster  chambers  may  limit  the  amount  of  achievable  thrust 
per  unit  area.  Other  factors  such  as  local  humidity  and  subtle  electrode  alignment  changes  greatly  affected 
the  maximum  attainable  ion  wind  velocity  [11]  indicating  that  EHD  propulsion  systems  may  not  be  very 
robust.  Obviously  many  other  factors  can  negatively  impact  the  performance  of  EHD-thrusters.  However, 
there  is  still  no  commonly  accepted  view  on  EHD-thruster,  e.g.,  whether  EHD-thruster  is  worthwhile  or 
not,  what  level  of  thrust  and  thrust  efficiency  can  be  obtained  from  EHD-thrusters,  etc.  Therefore,  the 
construction  of  an  ideal  model  of  EHD  thrusters  is  very  important  for  understanding  major  physical 
limitations  of  EHD  thrusters  in  "ideal  situations”. 
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In  this  paper  we  present  a  simple  one-dimensional  model  of  an  ideal  EHD  thruster.  A  description  of  the 
model  and  numerical  results  are  presented  in  Section  II  and  III,  respectively,  and  conclusions  are  given  in 
Section  IV. 


II.  Description  of  the  Model 

A  schematic  of  an  ideal  electrohydrodynamic  thruster  is  shown  in  Fig.  1 .  The  ions  are  “seeded”  at  the 
anode  and  are  absorbed  at  the  cathode.  Here  we  do  not  consider  the  mechanism  for  creating  the  ions.  They 
can  be  created,  for  example,  by  corona  discharge,  or  by  other  means.  The  neutral  gas  flow  is  coming  from 
the  anode  side  of  the  thruster.  Due  to  the  friction  forces  with  ions,  the  neutral  gas  flow  is  accelerated  in  the 
discharge  thruster  region  that  creates  the  thrust,  as  in  an  ordinary  air  turbine  jet. 

The  following  assumptions  are  made  in  the  model.  Fig.  1:  1)  the  width  of  the  thruster  is  much  larger 
than  the  length  of  the  thruster  L  ;  2)  all  parameters  of  incoming  neutral  gas  and  ion  flow  are  uniform  across 
the  thruster;  3)  the  anode  and  the  cathode  are  transparent  for  the  neutral  gas  flow;  4)  there  are  no  electrons 
in  the  discharge  region;  5)  all  parameters  of  the  thruster  are  independent  of  time  (i.e.,  we  assume  the  steady 
state  regime  of  thruster  operation);  6)  the  temperature  of  neutral  gas  flow  in  the  discharge  chamber  is 
constant.  The  sixth  assumption  relies  on  the  assumption  that  the  electrical  power  put  into  the  discharge  is 
relatively  small  so  that  the  heating  of  the  gas  flow  due  to  the  discharge  can  be  ignored.  It  follows  from 
assumptions  1-3  that  our  model  is  one -dimensional. 

The  equations  describing  one-dimensional  thrusters.  Fig.  1,  are 

dE 

sQ-  —  =  -e-ni  (1) 

dx 

j  =  e-nt  ■  (u  +  Vrfrijj  )  =  const  (2) 

na,Q  '  u0  =  na  ' 11  (3) 

du 

ma  '  na.O  ■U(y--  =  nre-E  (4) 

dx 

where  Eq.  (1)  is  Payson  law,  Eq.  (2)  is  electrical  charge  conservation  law,  Eq.  (3)  is  mass  conservation  law 
for  neutral  particles,  and  Eq.  (4)  is  momentum  conservation  law  for  ions,  and 

v drift  =  Mi  E  (5) 
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(6) 


e  ■  zin  M,- 
Mi  = - —  =  — 

mi  na 

In  Eq.  (2)  we  have  assumed  that  the  ion-neutral  collision  mean-free-path,  is  much  smaller  than  L  and, 
therefore,  the  ion  velocity  without  the  electric  field  has  to  be  about  equal  to  the  local  velocity  of  gas  flow.  It 
should  be  stressed  that  in  the  case  where  a  corona  discharge  is  used  to  ionize  the  incoming  gas  flow,  Eq.  (2) 
is  exact.  In  the  model  we  also  assume  that  M(  is  a  constant.  This  is  a  reasonable  assumption  as  long  as  the 
relative  velocity  between  the  gas  flow  and  the  ion  “wind”,  Vdn ^ ,  is  smaller  than  the  neutral  particle 


thermal  velocity,  Vja  ,  and  the  temperature  of  gas  flow  Ta  does  not  change  significantly  in  the  discharge 
region  of  the  EHD-thruster,  assumption  6.  It  should  be  stressed  that  in  the  case  of  strong  electric  fields, 
when  the  energy  that  an  ion  gains  from  the  electric  field  between  collisions  is  comparable  to  the 
temperature  of  the  neutral  gas  ( Vdlp  calculated  by  Eq.  (5)  becomes  larger  than  VTa ),  the  ion  drift  cannot 


be  described  by  Eq.  (5)  any  longer  [12],  It  is  worth  noting  that  this  situation  corresponds,  in  particular,  to 
the  case  of  the  electrical  breakdown  in  gasses.  The  system  of  Eqs.  (1)  -  (6)  is  essentially  the  system  of 
equation  of  Christensen-Moller  [7],  however  taking  into  account  dependence  of  ion  mobility  on  density  of 
neutral  gas. 

Integrating  Eq.  (4)  along  the  discharge  region  we  obtain  an  equation  for  the  thrust  per  unit  of  area  of 
the  discharge, 

L 

T  =  ma  ■  naQ  ■  Uq  ■  (ii]  —Uo)  =  e-j n;-  ■  E-dx  (7) 

o 


Multiplying  Eq.  (4)  by  u  and  then  integrating  along  the  discharge  region  we  obtain  an  equation  for  the  gain 
of  kinetic  energy  of  gas  flow  per  second  per  unit  of  area  of  the  discharge. 


Wk  =  ma  ■  na  Q  ■  u0 


M|  ll0 

T  ~2 


=  { iij  -e-E-u-dx 
0 


(8) 


The  electrical  power  deposited  in  the  discharge  per  unit  area  of  the  discharge  can  be  calculated  as 


Wj  =  j  ■  \ E  ■  dx  =  j  -U 
0 


(9) 


where 
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(10) 


U  =  \E-dx 
0 

The  system  of  Eqs.  (1)  -  (6)  can  be  reduced  to  the  following  system  of  two  ordinary  differential 
equations. 


dVdrift 

dx 


-  B  ■ 


VL, 


drif 


-A-- 


(V 'drift  +  u)  ■  U 1  (V ’drift  + 


^=B- 


V, 


drift 


dx  (V drift  +  U)  ■  U 


(ID 


(12) 


where  A  and  B  are  constants: 


A  =  - 


j '  M,- 


and  B  =  - 


J 


M,-  •  ma 


(13) 


^0  '  na, 0  '  w0 

We  have  numerically  solved  Eqs.  (1 1)  -  (13)  with  two  boundary  conditions: 

u(x  =  0)  =  tty  and  U  =Uq  (14) 

Now  let  us  consider  the  case  where  the  ion’s  drift  velocity  is  much  larger  then  the  velocity  of  neutral 
gas  flow, 

V drift  »  u  (!5) 

Dropping  u  in  Eq.  (2),  we  obtain  the  following  equation  for  E  , 

=  (16) 

dx  u  t  ■  E 

Assuming  further  that  //,  =  const  and  solving  Eq.  (16)  analytically  for  given  j  and  Eq  ,  we  obtain  the 
following  analytical  estimates  for  «/(x),  E(x) ,  U  ,  and  ymax : 

.  ^-1/2 


J 


e  ■  Mi  ■  Eq 


1-*. 

L 


V  J  max  J  J 


E  =  E, 


f  x  f  •  ^ 

i— — ■ 


V  imax  J ) 


1/2 


U  =  --Ea-L-dlM 

3  j 


1- 


/  .  \3/2^ 


1-- 


V  J  max  J 


(17) 


(18) 


(19) 
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where 


Jm  ax 


£0  '  Mi '  E()  _  9  ■  £q  ■  ft  ■  V 


2  L 


\-L3 


(20) 


It  follows  from  Eq.  (17)  that  the  current  cannot  exceed  the  maximum  space-charge-limited  current  jmax, 
the  Child -Langmuir  Law.  In  the  case  when  j  — »  0  we  obtain  that  U  =  Eq  ■  L  and  nt  =  0 ,  and  when 

2 

j  =  jmax  Eq.  (19)  yields  U  =  —  ■  Eq  ■  L .  Substituting  Eq.  (17)  -  (19)  into  Eqs.  (7)  and  (8)  and  setting  u  to 


Uq  in  the  right-hand-side  of  Eq.  (8)  we  obtain  analytical  estimates  for  the  thrust  and  the  thrust  efficiency. 


T 


jL 


Mi 


(21) 


W7/.  ii  ■  L  u 
~W~j  ~  ~Mi~^  ~ 


(22) 


Substituting  /max ,  Eq.  (20),  into  Eq.  (21)  we  obtain  an  estimate  for  the  maximum  thrust  that  can  be 


achieved  from  the  thruster. 


./max  ‘  ^ 

Mi 


9-Sq-U2 


■  L 


(23) 


As  one  can  see  from  Eq.  (23),  the  maximum  thrust  is  independent  of  ion  mobility  when  Vtlnf  »  u .  It 


follows  from  Eqs.  (22)  and  (23)  that  with  an  increase  in  the  thruster  voltage,  7max  increases  while  % 
decreases.  This  shows  that  achieving  both  targets,  high  thrust  and  high  thrust  efficiency,  at  once  can  be 
difficult. 


III.  Numerical  results 

The  exact  calculation  of  ion  drift  mobility  in  air  is  a  very  complicated  task;  notably  it  depends  on  ion 
composition,  humidity  of  air,  pressure,  and  other  factors.  Direct  measurements  of  average  natural  air  ion 

mobility  at  one  atmosphere  [13]  give  //,•  «  1 .39  •  1 0  4  m2/(V-sec)  with  a  standard  deviation  of  0.1-10  4 
m2/(V-sec).  Since,  the  ion  mobility  of  Ni  in  pure  nitrogen  gas  at  a  pressure  of  one  atmosphere  and  a 
temperature  of  300  K  is  1.37-10  m“/(V-sec)  [12]  and  close  to  the  measured  average  natural  ion  mobility 
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of  air,  we  consider  in  our  paper  the  case  of  pure  nitrogen  gas.  In  our  calculations  we  have  used 


M,-  =  1 .27  ■  1 024  [V/(sec-m)]  that  corresponds  to  drift  of  ions  N%  in  nitrogen  gas  at  gas  temperature  of 

94  n. 

300K.  Since  M,  is  weakly  dependent  on  the  temperature  (M,  =1.01  10“  m/(V-sec)  at  1000  K),  our 
assumption  that  M,-  =  const  is  well  satisfied. 

Figs.  2-13  present  the  result  obtained  by  numerically  solving  the  system  of  Eqs.  (11)  -  (14)  at 
maximum  current  (maximum  thrust).  Fig.  2  shows  maximal  current  densities  vs.  voltage  for  different 
velocities  of  neutral  gas  flow  Uq  ;  in  our  model  Uq  can  be  interpreted  as  the  speed  of  the  vehicle.  As  it  has 
been  expected  with  an  increase  in  applied  voltage,  ymax  increases  rapidly.  As  shown  in  Fig.  3  the 
maximum  thrust  of  the  EHD-thruster  increases  with  an  increase  in  the  voltage  (as  expected)  and  seems  not 
to  depend  greatly  on  the  speed  of  vehicle  as  in  the  case  of  »  u  ,  Eq.  (23).  Fig.  4  shows  that  thruster 

efficiency  increases  with  an  increase  in  the  speed  of  the  vehicle  and  decreases  with  an  increase  in  the 
applied  voltage;  these  trends  agree  with  an  estimate  of  %  for  the  case  of  V(/lip  »  u ,  Eq.  (22)  and 

qualitatively  agree  with  the  Bondar-Bastein  model  [9], 

To  investigate  how  the  length  of  EHD-thruster  affects  the  thruster  performance  we  have  calculated 
parameters  of  the  EHD  thruster  for  Uq  =  50  m/s  and  L  =  1  m  for  different  /max ,  Figs.  5  -  10.  As  shown  in 
Fig.  5  the  discharge  voltage  increases  very  non-linearly  with  the  length  of  the  capillary  for  fixed  magnitude 
of  ymax.This  behavior  of  U  vs.  L  can  be  explained  by  significant  non-linear  dependencies  of  the  electric 
field  on  discharge  for  currents  close  to  y'max ;  see  an  analytical  expression  for  U  in  Eq.  (20)  for  the  case  of 
V drift  »  u  ■  With  an  increase  in  L  the  thrust  increases  more  or  less  linearly.  Fig.  6,  but,  the  thrust 
efficiency  decreases  sharply.  Fig.  7.  We  have  also  calculated  AV  ,  one  of  the  important  parameters  of 
thruster  performance.  Fig.  8.  As  one  can  see  AV  reaches  about  4.5  m/s  for  L  =  100cm  and  /max  =  0.05 
A/m2.  Since  the  ion  current  cannot  exceed  the  space-charge-limited  current  jm ax ,  the  ion  density  in  the 
discharge  is  very  small.  Fig.  9. 

To  investigate  the  performance  of  the  thruster  at  different  altitudes  we  have  also  calculated  the  thrust 
and  thruster  efficiency  for  H  =  20  km.  As  shown  in  Figs.  10  and  11,  both  thrust  and  thrust  efficiency 
decrease  with  an  increase  in  H  .  The  reason  is  that  ion  mobility  increases  with  a  decrease  in  gas  density. 
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—1  —1/2 

Eq.  (6),  leading  to  a  decrease  in  7'max  and  %  ( 7max  cc  /i-  and  x  oc  /t,  ,  Eqs.  (22)  and  (23)).  Thus, 

since  the  drag  and  the  thrust  are  approximately  proportional  to  gas  pressure,  using  EHD-thrusters  at  high 
altitudes  does  not  make  much  sense. 

To  investigate  to  what  extent  the  simple  analytical  Eqs.  (22)  and  (23)  are  applicable,  we  have 
calculated  the  ratios  of  analytically  predicted  thrust  and  thruster  efficiency  to  numerically  calculated  rmax 
and  x  >  Figs.  12  and  13  respectively.  Fig.  13  also  shows  the  Bondar-Bastein  estimated  thruster  efficiency 
[9], 


X  = 


«o 


«0  + 


Mi  'U 


(24) 


which  coincides  with  Eq.  (22)  when  V^ri^  »  u .  As  one  can  see  Eqs.  (23)  and  (24)  give  reasonable 
estimates  for  the  7max  and  x  in  almost  all  regions  of  /nlax ,  while  Eq.  (22)  works  only  for  large  /max 
corresponding  to  large  Vdrift . 

The  increase  in  temperature  of  gas  flow  due  to  Ohmic  heating  by  the  discharge  can  be  estimated  as 


2-(Wj-Wk-nrvjift) 
3 -k-na  ■  HQ 


(25) 


As  shown  in  Fig.  14  ,  A Ta  increases  with  a  decrease  in  speed  of  the  vehicle  and  with  an  increase  in  the 
length  of  the  thruster;  this  makes  perfect  sense.  The  faster  neutral  gas  streams  through  the  EHD-discharge 
chamber,  the  smaller  A Ta  gets.  Since  A Ta  is  quite  small  even  for  L  =  lm,  the  sixth  assumption  in  Section 
II  is  well  satisfied. 

Table.  1  shows  the  ratio  of  the  ion  drift  velocity  at  the  anode  (where  electrical  field  is  largest)  to  the 
thermal  velocity  of  neutral  gas  flow  and  U  l(P  ■  L )  for  different  altitudes,  gas  flow  velocities,  and  lengths 
of  the  EHD-thruster.  As  was  mentioned  in  Section  II,  when  t//(.(/,  /VTa  approaches  unity,  Eq.  (5)  for  ion 

mobility  becomes  invalid.  Since  the  electrical  breakdown  voltage  of  nitrogen  gas  is  about  1.4- 106 
V/(m- Atm)  [12]  (and  is  close  to  the  electrical  breakdown  of  air  [14]),  the  electrical  breakdown  most  likely 
appears  before  Eq.  (5)  became  invalid.  Table  1.  Therefore,  the  obtained  results  in  the  case  when 
Vdrift  /V Ta  >  1  should  be  considered  from  an  illustrative  point  of  view  only. 
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IV.  Conclusions 


We  presented  a  one -dimensional  model  of  an  ideal  EHD-thruster.  We  showed  that  the  ion  current  in  an 
ideal  EHD-thruster  has  a  principal  limit;  it  cannot  exceed  the  space-charge-limited  current  ymax,  and 
therefore,  the  thrust  of  such  a  thruster  cannot  exceed  the  thrust  corresponding  to  this  current.  We  also  show 
that  thrust  and  thrust  efficiency  “work”  against  each  other  in  EHD -thrusters:  the  larger  the  thruster  current 
(larger  thrust)  the  smaller  the  thruster  efficiency.  We  also  show  that  at  high  altitudes  the  performance  of 
EHD-thrusters  (thrust  and  thrust  efficiency)  drops  very  fast,  and  therefore,  using  EHD-thrusters  at  altitudes 
larger  than  5  km  is  apparently  unrealistic.  The  model  shows  that  maximum  thrust  cannot  exceed  20  -  30 
N/m2  at  sea  level  even  at  breakdown  voltage;  thus,  this  also  does  not  look  very  encouraging.  Since  our 
model  is  an  ideal  and  does  not  take  into  account  numerous  other  factors  that  can  negatively  impact  the 
performance  of  EHD-thrusters  even  more,  we  may  conclude  that  the  idea  of  using  the  Electrohydrodynamic 
ionic  wind  pump  effect  for  thruster  application  is  rather  dim. 
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Fig.  1.  Schematic  (not  to  scale)  of  electrohydrodynamic  thruster. 


Ideal  EHD  Thruster,  N2,  H  =  0  km,  L  =  10  cm 


Ideal  EHD  Thruster,  N2,  H  =  0  km,  L  =  10  cm 


U  [Volts] 

Fig.  2.  Maximum  current;  solid  line:  Uq  =  3  m/s, 
dashed  line:  Mq  =  m/s,  dotted  line:  Uq  =  50  m/s. 
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Fig.  3.  Thrust  at  maximum  current;  solid  line:  Uq  =  3  m/s, 
dashed  line:  Uq  =  10  m/s,  dotted  line:  Mq  =  50  m/s. 
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Ideal  EHD  Thruster,  N2,  H  =  0  km,  L  =  10  cm 


U  [Volts] 

Fig.  4.  Thrust  efficiency  at  maximum  current;  solid  line: 

Uq  =  3  m/s,  dashed  line:  u g  =  10  m/s,  dotted  line:  u q  =  50  m/s. 


Ideal  EHD  Thruster,  N2,  H  =  0  km,  V0  =  50  m/s 


jmax  [A/m  ] 


Fig.  5.  Voltage  at  maximum  current;  solid  line: 
L  =  10  cm,  dashed  line:  L  =  100  cm/s. 
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Ideal  EHD  Thruster,  N2,  H  =  0  km,  VO  =  50  m/s 


jmax  [A/m  ] 

Fig.  6.  Thrust  at  maximum  current;  solid  line: 
L  =  10  cm,  dashed  line:  L  =  100  cm/s. 


Ideal  EHD  Thruster,  N2,  H  =  0  km,  V0  =  50  m/s 


jmax  [A/m  ] 


Fig.  8.  Delta-V  at  maximum  current;  solid  line: 
L  =  10  cm,  dashed  line:  L  =  100  cm/s. 


Ideal  EHD  Thruster,  N2,  L  =  10  cm,  V0  =  50  m/s 


Jmax  [A/m2] 


Fig.  10.  Thrust  at  maximum  current;  solid  line: 
H  =  0  km,  dashed  line:  H  =  20  km. 


Ideal  EHD  Thruster,  N2,  H  =  0  km,  VO  =  50  m/s 


jmax  [A/m  ] 

Fig.  7.  Thrust  efficiency  at  maximum  current;  solid 
line:  L  =  10  cm,  dashed  line:  L  =  100  cm/s. 


Ideal  EHD  Thruster,  N2,  H  =  0  km,  V0  =  50  m/s 


jmax  [A/m  ] 


Fig.  9.  Ion  density  at  maximum  current;  solid  line: 
L  =  10  cm,  dashed  line:  L  =  100  cm/s. 


Ideal  EHD  Thruster,  N2,  L  =  10  cm,  V0  =  50  m/s 


jmax  [A/m  ] 


Fig.  11.  Thrust  efficiency  at  maximum  current;  solid 
line:  H  =  0  km,  dashed  line:  H  =  20  km. 
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max, estimate  '  1  max, model 


Estimates  vs.  Model: 

H  =  Okm,  L  =  10cm,  u0  =  50m/s 
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Fig.  12.  Ratio  of  the  estimated  thrust  to  the  numerically 
calculated  thrust  at  maximum  current. 


Ideal  EHD  Thruster,  N2,  H  =  0  km 


jmax  [A/m2] 


Fig.  14.  An  estimation  for  A Ta  ;  solid  lines:  L  =  lm, 
dashed  lines  L  =  10cm 


Estimates  vs.  Model: 
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Fig.  13.  Ratio  of  the  estimated  thrust  efficiency  to  the 
numerically  calculated  thrust  at  maximum  current; 
solid  line:  Eq  (22),  dashed  line:  Eq.  (24). 
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Table  1.  Vdrip  /VTa  and  E/P  [Volts/(m- Atm)| 

H  =  0 km  L  =  0.  bn  L  =  bn  H  =  20 hn  L  =  0.  bn 

Vq  =  3m  Is  V0  =  1  Oin  /  s  V0  =  50m  /  s  V0  =  50m  /  s  V0  =  50m/  s 


A/ Hi2 

V drift  ^ Ta 

E/P 

V drift  ^ Ta 

E/P 

V drift  ^ Ta 

E/P 

V drift  ^ Ta 

E/P 

£ 

E/P 

0.001 

0.07 

8.63xl04 

0.12 

1.80x10s 

0.13 

2.09x10s 

0.34 

5.02x10s 

0.61 

7.85x10s 

0.005 

0.022 

3.10x10s 

0.29 

4.56x10s 

0.29 

4.58x10s 

0.89 

1.36xl06 

1.53 

2.04xl06 

0.01 

0.34 

5.02x10s 

0.42 

6.06x10s 

0.40 

6.15x10s 

1.29 

2.02x1 06 

2.22 

2.98xl06 

0.05 

0.89 

1.36xl06 

0.92 

1.45xl06 

0.74 

1.09xl06 

2.97 

4.70x1 06 

5.12 

6.99xlOb 
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